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The first order QED corrections to the polarized muon decay spectrum are consid- 
er ' ered. The exact dependence on electron and muon masses is kept. Numerical results 
■ are presented. 
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1 Introduction 



Since tlie discovery of muon in 1936, experimental and theoretical investiga- 
^ ■ tions of its properties became an important part of the elementary particle 

physics. The very accurate measurements of the properties provide serious 
checks of the Standard Model and give a possibility to look for new physics 
in low-energy experiments. Because of the more and more precise experi- 
mental facilities and techniques, calculations of radiative corrections become 
unavoidable. We need them to obtain theoretical predictions with the required 
accuracy. 

In this paper we will consider radiative corrections to the polarized muon decay 
spectrum. This work was initialized by the experiment TWIST [1,2], which 
is going to measure the spectrum with the error level of the order of 10~^. 
We will speak here only about the first order quantum electrodynamic (QED) 
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correction and concentrate on the effect of the non-zero ratio of the electron 
and the muon masses. Higher order corrections to the muon decay spectrum 
wiU be considered elsewhere [3,4]. A review about the muon properties and 
possible non-standard effects in the muon decay can be found in Refs. [5,6]. 



2 The Born level 



At the Born level within the Fermi Model the differential width is described 
by the formula [7,8]: 



j2pBorn C'^ 

/^°-(x)=3-2x + |(3x-4)(l-/3^). 



-^1""^' Ee = ^x, c = cose, 



where nie and are the electron and the muon masses; 9 is the angle between 
the electron momentum and the muon polarization vector (c — > — c for the 
decay); ^ is the degree of the muon polarization; /3 is the electron velocity 
in the muon rest reference frame; is the electron energy; x is the electron 
energy fraction, 

ml 

In the massless limit {m^ 0) we have 

/Born^^) ^ = 3 - 2x, g^'^^^ix) ^ go{x) = 1 - 2x. (3) 
The integration over the energy fraction gives 

Jdx x'f3f°^^ix) = -F(p), Jdx x'fJg^'^^-ix) = -G{p), 

F{p) = l-8p- 12p2 In p + 8p2 - p^ 

G{p) = -i + fp'/^ - 30p^ + 32pV2 _ + ^4 
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Function F{p) is relevant for the total decay width at the Born level: 



r^°- = jdcTo^^roFip). (5) 
-1 

Function G{p) contributes to the forward-backward asymmetry of the decay: 
= / dc ce - / dc ce Fo^ = eroG'(p). (6) 

-1 

Within the Standard Model the muon decay happens due to weak interaction 
of leptons and li^-bosons. The Fermi Model corresponds to the limiting case 
of the infinite 14^-boson mass. We follow here the definition of the Fermi cou- 
pling constant Gp as discussed in Ref. [9]. That means, all weak effects are 
incorporated into the coupling constant, and QED radiative corrections have 
to be calculated within the Fermi Model. I accept and support this approach. 
Originally in the literature, the constant is defined in a different way [10-12], 
so that the first order effect in the muon and the H^-boson mass ratio gives 



m. 



(7) 



But in any case, in studies of the muon decay spectrum, it is natural to 
use the constant directly defined from very precise experiments on the muon 
lifetime. We have to note here, that although the Fermi lagrangian itself is 
not renormalizable, QED corrections to the process under consideration can 
be shown to be finite [13] at all orders of the perturbation theory. 



3 The Exact First Order QED Corrections 



Here we will consider the O (a) QED correction to the muon decay spectrum 
with keeping the exact dependence on the electron and the muon masses. The 
result is obtained by means of the standard technique. The contributions from 
virtual, soft, and hard photons were evaluated separately: 

d^r*^^^ d^r^'""* (;j2-pSoft ^2-pHard 

dxdc dxdc ~^ dxdc ~^ dxdc ^ ^ 
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To be short, we give here only the simple formula for the soft photon contri- 
bution: 



j2FSoft j2FBorn 

— = — ^Soft^ /gx 

dxdc dxdc ' 



2}, h = ^^Y^: (10) 



where Ae is the maximal energy of a soft photon (A -C 1); A is a fictitious 
photon mass; functions Li2 (x) and ({n) are defined in the Appendix. 

The auxiliary parameters A and Ae cancel out in the total sum of the three 
contributions: 



^ = To x'P^ (Mx) + cCgiix)) , (11) 



+ 



p 2xP 
4 - x2(l - 2 - x{l - P) 



m. 



xf3 



1/. l + (3\(5x^ 



+ ^(^^i^ + 21n. + 21n^j{|^(l-,^)3-^(l-,^)^ 
+ ^(3 - 12/3 + /3^)(1 - /3^) + + 2/3 + (1 - + 



32 



+ ^[-20 - 12/3 - 19(1 - /32)] + ^ + —\ 



X 



j 



+ (^In X + In 

I yBorn 



1+/3 
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x\l- P^) + 2xiP^ -3) + 3 



[X) 

22 



2 

11 22 2 2 



-— x(l - /3') + — /3' - 
18 ^ ^ ' 27 9 



-(-i^*4-^)-f<-''vi, 



(12) 



ii,„4_,,,+i„i±^+i„?-(i-« 



mf 



+ 2 In 



/5 



+ 4 In 



2/3 

2 -a;(l -/3) 
2^ 



2/3 
In X + In 



+ 



In — - — 2 In X 



1 + /3 



4 



,(l + /3)(2-a;(l -/?)), 



(13) 



2 



+ a; 



2 2 



x(^x^(l-/?^)-4x + l)+i,(ln-^-21n 



X 



+ a; 



a; 

'48 



.3^3 _ i 
2^ 4^ 4 



4/3+(l-/3^)(-^/3^-4/3-^) 



.f-2. 



+ (l-/32)(2/3-2) + -(-6 + (l-/3^)) + 



a; /17 

+ 



2\2 



144/^2 

- 195/?) 



,x2(l-/5^) 



; 324 V /3 

2 _2 



1 /595 10 
— - 1923/? + 



10 55 203 

,t(1 - /3 ) - 77 + 



54 162/32, 
1 



L = ln— |, ^ ml(l - x) + ml. 



J 3x(3 a;2/3' 



(14) 



Again the integration over the energy fraction and the angle gives us correc- 
tions to the decay width and the forward-backward asymmetry: 



(15) 



i^i(p) = (1 - P^) (f - + f P^) - pin p (20 + 90p - 

+ ^p^) - In^ p(36 + p^) - (1 - /) (y - f P + y P^) 

X ln(l - p) + 4(1 + 30p2 + p^) lnpln(l - p) + 6(1 + 16p2 + p^) 
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X [Li2 (p) - C(2)] + Up^'\l + p) ( 3C(2) - 2Li2 (Vp) 



+ 2Li2 (-^) - In pin , (16) 



^/x 472, , x/1271 47 2959 

4657 2 ,,5/2 21 ,\ , ,/100 52 

268 o 20 A , . „2048 



+ ^P^ + fp^)+ln(l-Vp)(l-VP) 

,/ 2035 2048 ^ 1987 512 o,. 
+ ln(l - VP) (1 - P) (-- ^ + ^ VP - + —P''" 

637 2 17 3\ , ,^ ^,608 , ,^ ,/608 
" "9"^ + j " ^"^^^ " Vp)-^ + lnp(l - P)(— 

614 ^ 623 3/2 967 2 ^5/2 1 3 

- Vp + -^P - 56p'/' + — p' + 2p5/2 + -p3 

/2 S2 

+ In p[ln(l - Vp) + 2 ln(l + Vp)] ( 3 + 8p^ - y p^ + 2p^ 

/14 
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+ p2 In^ p(7 + 16p - -p2) + (C(2) + 2Li2 (-^p)) (y + 16p 

+ 32p2 - yp3 + 6p^) - p2^ In^ p + 16p2 lnp[2C(2) - 2U2 {-^/p) 

- Li2 (Vp)] + 16p2[2Li3 (Vp) + 12Li3 (-^P) + 7C(3)]. (17) 

The definitions of the special functions, used above, are given in the Appendix. 

The expression for -Fi(p) coincides with the one, received in Ref. [14], starting 
from a differential distribution of a different kind. We reproduced here the 
formula for the sake of completeness and for a comparison with Gi(p). Both 
the functions Fi{p) and G'i(p) are vanishing in the limit p ^ 1, because of the 
vanishing phase space volume. 

It is interesting to note, that Gi(p) contains terms of the first power in the 
electron and the muon mass ratio, while Fi{p) does not (see the discussion 
about the odd mass terms in Ref. [9]). The linear mass terms in Gi(p) can be 
clearly seen from the expansion ^ of the exact formula: 

f ^^^^ " + K"^^ + ^^^^^^ + ^) 

+ P'^' + ^ 1- ^) + P" + 32C(2) + 112C(3) - 16 In p 



The expansion of Fi{p) can be found in Ref. [14]. 
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oo^/oM ^,1 2 2 ,o \ ./2/ 95624 1232, ' 
+ 32C(2) Inp + 7 InV - 3 ln=^ pj + p'^' ("^^ ^^^""^ 

o/5662 16^,^, 698, ,^,2 
+ P [-^ - yC(2) - — lnp+ 161nV 

7/2/ 134248 512, \ 4, 

One can check that the term with |me/m^|^ = ^/p is coming from the integra- 
tion over the region of large (close to the upper limit) values of the electron 
energy fraction. The absolute value appears from the integration, it can be 
important for analytical continuations. Analogous terms with the first order 
mass ratio have been found in the forward-backward asymmetry in the process 
of electron-positron annihilation into heavy leptons [15]. 

The incorporation the first order correction gives the following formula for the 
spectrum: 



^2prad.corr. ^ ^2pBorn ^ ^2p(l) ^ ^ 

dxdc dxdc dxdc 



4 Comparisons and numerical results 



The agreement with the massless formulae [16] for /i(x) and gi{x) is checked. 
Our calculations for the massive case agree in part with the known results. 
Namely, function f\{x) coincides with the result of Ref. [17], where a mistake 
should be corrected according to Ref. [16] (see also Appendix C in Ref. [9]). 
The same function does agree with the one given in Ref. [18], while the the 
integral over the real photon phase volume has not been taken analytically 
there. Moreover, a confirmation of the given formula for fi{x) comes from 
the comparison of the result of its integration over the energy fraction with 
Ref. [14], as mentioned above. 

Function gi{x), which describes the polarized part of the decay spectrum, 
could have been calculated long time ago on the same basis as fi{x). Neverthe- 
less, our results for gi{x) and Gi(p) are new. A partial comparison between the 
virtual (loop) diagram contributions shows an agreement with the correspond- 
ing quantity from the calculation of the top-quark decay spectrum [19] Q. 

In Table 1 we show the effect of the mass terms at the lowest (Born) level 
and at the level of the first order correction for different points of the energy 
spectrum; 

^ Because of different choices of observables, the comparison of the total spectra is 
impossible. 
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Table 1 

The effect of tlie finite electron-muon mass ratio versus tlic electron energy fraction. 



X 


ho{x) 




jfBorn 
"m 




hi{x) 








c 


= 0, 








0.05 


0.00725 


0.00711 


-194.5 


4.11481 


4.10454 


-16.45 


0.1 


0.02800 


0.02786 


-49.52 


5.95508 


5.95444 


-0.266 


0.2 


0.10400 


0.10387 


-12.80 


8.68399 


8.68517 


0.132 


0.3 


0.21600 


0.21587 


-5.796 


10.3054 


10.3067 


0.071 


0.5 


0.50000 


0.49989 


-2.105 


8.66761 


8.66871 


0.026 


0.7 


0.78400 


0.78391 


-1.088 


-1.55489 


-1.55410 


0.012 


0.9 


0.97200 


0.97193 


-0.742 


-25.6678 


-25.6670 


0.010 


0.99 


0.99970 


0.99963 


-0.702 


-67.6027 


-67.6011 


0.019 


0.999 


1.00000 


0.99993 


-0.702 


-107.665 


-107.663 


0.026 








= 1. 


C = 1 






0.05 


0.00950 


0.00928 


-236.8 


3.6880 


3.6657 


-27.32 


0.1 


0.03600 


0.03579 


-59.02 


5.2896 


5.2850 


-1.490 


0.2 


0.12800 


0.12781 


-14.51 


7.4177 


7.4179 


0.027 


0.3 


0.25200 


0.25184 


-6.186 


7.8913 


7.8925 


0.054 


0.5 


0.50000 


0.49991 


-1.871 


2.7579 


2.7589 


0.024 


0.7 


0.58800 


0.58796 


-0.659 


-8.1130 


-8.1129 


0.001 


0.9 


0.32400 


0.32400 


-0.150 


-13.120 


-13.121 


-0.022 


0.99 


0.03920 


0.03920 


-0.013 


-3.2805 


-3.2807 


-0.053 


0.999 


0.00399 


0.00399 


-0.001 


-0.4949 


-0.4949 


-0.081 



^Born ^ -^Q^ 



- 1 



^Born 

ho{x) 

ho,i{x) = fo,i{x) + c^go,i{x), 



= 10^ 



a hi{x)-K^^\x) 



271 ho{x) 

where the functions and g^^^^(x) can be received by applying 

rrie = in the general expressions for fi{x) and gi{x) everywhere, except 
the argument of the large logarithm L. In practice, the massless functions 
can be taken from Ref. [16]. The argument of the functions with the exact 
mass dependence is rescaled: x — x ■ Xmax- One can see, that the effect of the 
mass terms in the 0{a) order is below the 10~^ precision tag (sec S^^^) for 
experimentally [1,2] preferable values of the electron energy fraction x > 0.3. 



In Table 2 we show the effect of the mass terms in the integrated quantities: 
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Table 2 

The effect of the finite mass ratio in the total decay width and in the forward- 
backward asymmetry. 



p 


m 


np) 


rBorn 
Op 


Fi(0) 


Flip) 




^ = 2.34 • 10'^ 
^ = 3.54 • 10-3 


1.0000 
1.0000 


0.9998 
0.9726 


-1.871 
-274.3 


-3.6196 
-3.6196 


-3.6152 
-3.3367 


0.051 
3.286 


P 


G(0) 


G{p) 


^Born 


Gi(0) 


Gi{p) 


.(1) 


^ = 2.34 • 10-5 
^ = 3.54 • 10-3 


-0.3333 
-0.3333 


-0.3333 
-0.3314 


-0.036 
-56.71 


1.9634 
1.9634 


1.9502 
1.7651 


0.460 
6.908 



rBorn _ 1 n4 Y EiEl _ 1^ 
'^^ VF(0) V' 

rBorn _ I . (9!^^ 1^ 



Op — 



ill 

10 



^ Fi(p) - Fi(0) 
27r F(0) 
4 a G'i(p)-Gi(0) 



27r 



(21) 



The effect due to finite electron mass in function is approaching the 

10"^ level (for muon decay). 



5 Conclusions 



The present calculations can be easily extended for the case, where the final 
electron polarization is measured. But, for the moment, the experimental pre- 
cision there does not call for an account of small mass terms in the theoretical 
predictions. 

In this way, we considered the first order QED radiative correction to the 
muon decay spectrum. Formulae for the polarized part of the spectrum {gi{x) 
and Gi) with the exact dependence on the electron mass are received for the 
first time. The results presented are relevant for modern and future precise 
measurements of the muon decay spectrum. Our formulae are valid for the 
leptonic r decays as well. 
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Appendix 

Definition of special functions 



The Riemann ^-function and the pohlogarithm functions are defined as usu- 
ally: 



°° 1 vr^ 
Ci^) = E ' C(2) = - , C(3) = 1.20205690315959 . . . 

k=i " 

Lb(,) = -/dx!!^<l^, Li3M=/d.HiM. (A.1) 

J X J X 



The following identities can help in numerical evaluations and estimates: 



Li2 (1) = C(2), Li2 (-1) = -^C(2), Lis (1) = C(3), 
Lis (-1) = -^C(3), Li2 (y') = 2(Li2 (y) + Li^ (-y)), 



4 

Ligfi/^) =4(Li3(i/) + Li3(-y)). (A.2) 
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